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VIDEO CODING METHOD OF EXPLOITING THE TEMPORAL REDUNDANCY BETWEEN SUCCESSIVE 
FRAMES 



5 The present invention relates to a video coding method of exploiting the temporal 
redundancy between successive frames in a video sequence. 



Efficiently encoding video or moving pictures relies heavily on exploiting the 
temporal redundancy between successive frames in a sequence. Based on the as- 

10 sumption that local motions are slow with respect to the temporal sampling pe- 
riod, several techniques have been proposed for efficiently removing this redun- 
dancy. The most successful and acclaimed method is block-based motion predic- 
tion, heavily used in nowadays standards such as MPEG4 and H.26L. Roughly 
speaking, these compression scheme predict a frame in the sequence based on the 

15 knowledge of previous frames. The current frame (or predicted) is cut into blocks 
of fixed size and the best matching block is searched in the reference frame. Dis- 
placement vectors are then encoded so that the decoder can reconstruct the predic- 
tion of the current frame from the previously decoded frame(s). As the block- 
based prediction is not accurate enough to encode perfectly the current frame, the 

20 error between the original and the predicted frame is encoded separately. This is in 
general referred to as texture coding or motion residual coding. The main draw- 
back of this method lies in the blocky nature of the prediction mechanism, which 
gives rise to very noticeable blocky artefacts at low bit rates. Moreover such a 
system, while well suited for wide translational motions, is unable to cope with 

25 locally complex movements or even global geometric transformations such as 
zoom or rotation. Finally, block based motion prediction is not able to follow 
natural features of images since it is stuck in a fixed framework based on artificial 
image primitives (blocks). 
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This invention propose to solve the above-mentioned problems by introducing a 
new paradigm for dealing with spatio-temporal redundancy. 

The method according the invention is defined in claim 1 . In the dependent claims 
5 various embodiments are proposed. 

The invention will be described with the help of accompanying representations. 

Figure 1: shows the progressive reconstruction of an I frame with 50, 100, 150 and 
10 200 atoms. 

Figure 2 : shows a flow-chart of the I- frames codec, where the atoms are used to 
transmit the frame. 

Figure 3 : shows a flow-chart of the I-frames codec, where the atoms are estimated 
15 from the coded image both at the encoder and decode. 

Figure 4 : shows three successive schematic updates of basis functions inside a 
sequence of frames 

20 Figure 5 : shows a flow-chart of the P-frames coder, using atom based motion 
prediction. 

Figure 6 and 7: show the encoding of the Foreman sequence: I- frame, predicted 
frames #20, 50 and 99. 

25 

Figure 8: represents PSNR along the encoded Foreman sequence as a function of 
the number of predicted frames. 
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The proposed method is composed of two main parts. In a first step, a geometric 
model of a reference frame is built. In the second step, the model is updated by 
deformations in order to match successive frames. 

5 A reference frame (or intra frame, I-frame) is first decomposed into a linear com- 
bination of basis functions (atoms) selected in a redundant, structured library (see . 
P. Vandergheynst and P. Frossard, Efficient image representation by anisotropic 
refinement in matching pursuit, in Proceedings of IEEE ICASSP, Salt Lake City 
UT, May 2001, vol. 3, the content of which is incorporated herein by reference) 

0 

n=0 



In this equation, I(x f y) is the Intensity of the I-frame represented as a function giv- 
ing the gray level value of the pixel at position (x,y). c n are weighting coefficients 
15 and g (x, y) are the atoms involved in the decomposition. These atoms are image 

patches generated by a simple mathematical formula expressing the gray level 
value at pixel position (x,y). The formula is built by applying geometric transfor- 
mations to a function g(x,y) that we call a generating mother function. The pa- 
rameters of these transformations are concatenated in the vector of parameters y. 
20 Examples of possible transformations are translations, rotations or dilations. They 
act on the generating mother function by change of variable, for example: 



Translations: gbfcy) = g(x-b Jt y-b2) 

Dilation: gafcy) = a 1 g(x/a,y/a) 

25 Translations, anisotropic dilations and rotations: 
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cos £ (x - b } ) - sin & (y - b 2 ) 



a, 



g r ( x >y) = g( x n>y«)> where 



sind( 



x-b l ) + cosS(y-b 2 ) ' 
a 2 



y n = 



and y=[ai,a2,bi,b2,d], are the parameters of this transformation. 



Generating mother functions are chosen almost arbitrarily. Their properties can be 
5 adapted to the specific application. A possible example is to select an oscillating 
function of the form: 



The decomposition can for example be accomplished using a Matching Pursuit 
10 (see S. Mallat and Z. Zhang, Matching Pursuits with time- frequency dictionaries, 
IEEE Transactions on Signal Processing, 41(12):3397-3415, December 1993, the 
content of which is incorporated herein by reference). Matching Pursuit (MP) is a 
greedy algorithm that iteratively decomposes the image using the following 
scheme. First the atom that best matches the image is searched by maximizing the 



15 scalar product between the image and the dictionary (l|g r ), and a residual image 

is computed: 



g(x,y) = (l-x 2 )exp 




2 \ 



J 




20 Then the same process is applied to the residual: 



R,={R,|g r X+R 



and iteratively: 
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R„K R n|gr 0 )gr 0 +R n + i- 

Finally this yields a decomposition of the image in terms of a sum of atoms: 

5 

I W = Z(R n |gr B )gJ X ) + R NW. 

The basis functions (atoms) are indexed by a string of parameters y n representing 
geometric transformations applied to a generating mother function g(x,y). This 
10 index can be seen as a point on a manifold. The set of geometric transformations 
is designed in such a way that the total collection of basis functions (atoms) is a 
dense subspace of L 2 (R 2 ) , i.e. any image can be exactly represented. 

This part of the method expresses the I-frame as a collection of atoms that can be 
15 seen as geometric features such as edges or parts of objects which are very notice- 
able by the human eye. These basic primitives hereafter referred to as atoms, form 
a primal sketch of the image. The atoms are modelled and fully represented by the 
set of coefficients and parameters {c n , y n9 n = 0,...,N-l} , where c„ is the coeffi- 
cient and y n is a vector of parameters. 

20 

There are two ways to handle the coding of the I-frames. The first one is to esti- 
mate the atoms of the original frame. The atoms modelling the I frame are then 
quantized, entropy coded and sent in the bitstream. The process of quantizing an 
atom corresponds to the quantization of its coefficient and parameters. The atoms 
25 are also stored in memory in order to be used for the prediction of the next frames. 
A flowchart of this procedure is shown in Figure 2, where the atoms are used to 
transmit the frame. Dotted lines represent the coefficients and parameters of the 
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atoms, the bold solid line represents a video frame and the light solid line repre- 
sents bitstreams. The flow-chart should be understood as follow: The input frame 
is decomposed into atom using the MP algorithm. The atoms are then quantized 
and entropy. They are also de-quantized and stored in memory for the prediction 
5 of the next frames. Those de-quantized atom parameters are used to reconstruct 
the image, with the help of the generating mother function (which is known both 
by the encoder and decoder). Finally the difference between the original frame and 
the reconstructed one is computed and encoder using the frame coder. 

10 Note that the figure includes an optional step, that encodes the motion residuals 
(or texture), which is the difference between the original frame and the one recon- 
structed using the atoms. This encoding can be used to further increase the quality 
of the decoded image up to a lossless reconstruction. 



15 The second way of handling the I- frames is more conventional. The original frame 
is encoded and transmitted using any frame codec. Then the atoms are estimated 
from the reconstructed frame, both at the encoder and at the decoder. Finally those 
atoms are stored in memory for the prediction of future frames. The flowchart of 
this procedure is shown in Figure 3, where the atoms are estimated from the coded 

20 image both at the encoder and decoder. Dotted lines represent the coefficients and 
parameters of the atoms, the bold solid line represents a video frame and the light 
solid line represents bitstreams. The flow-chart should be understood as follow: 
The input frame is first encoder with the frame coder and send to the decoder in a 
bitstream. It is also decoded, in order to get the exact same frame that will be 

25 available at the decoder. This frame is decomposed into atoms by the MP algo- 
rithm. Finally those atoms are stored in memory to be used for the prediction of 
the next frames. 
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The second step of the method consists in updating the image model (the set of all 
atoms) in order to take into account the geometric deformations that have occurred 
between the reference and the current frame. Clearly, since the model is based on 
geometric transformations, updating its atoms allows for adapting to smooth local 

5 distortions (translations, rotations, scales are common examples). In order to com- 
pute this update, we assume that the deformed model is close enough to the refer- 
ence model We thus have to search for new atoms parameters in the proximity of 
the previous solution. This is performed by means of a local optimization proce- 
dure trying to minimize the mean square error between the updated model and the 

10 current frame (Figure 4, where three successive schematic updates of basis func- 
tions (atoms) inside a sequence of frames are represented). The updated atom pa- 
rameters are then the solution of: 



ArgMinj E^g,.- 1 * 

{c 0 .r n > HI n 

15 

where the optimization method for frame I t at time t is initialized with the atom 
parameters corresponding to the solution at time t - 1 or to the reference frame (the 
I frame) in order to avoid error propagation. This problem is a non-convex, non- 
linear, differentiable optimization problem (see Dimitri P. Bertsekas (1999) 

20 Nonlinear Programming: 2nd Edition. Athena Scientific 
http://www.athenasc.com/nonlinbook.html, , the content of which is incorporated 
herein by reference), which can be solved using various algorithms such as quasi- 
Newton methods, combined with line-search or trust-region globalization tech- 
niques (see Conn A., Gould N. & Toint Ph. (2000) Trust Region Methods. SIAM. 

25 http://www.fundp.ac.be/~phtoint/pht/trbook.html, the content of which is incorpo- 
rated herein by reference ), in order to identify a local optimum. 
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The difference between the original atom parameters and the updated ones is then 
computed and sent together with updated atom coefficients. Quantization and en- 
tropy coding can then be performed (P. Frossard, P. Vandergheynst and R.M, Fi- 
gueras y Ventura, Redundancy driven a posteriori matching pursuit quantization. 

5 ITS Technical Report 2002, the content of which is incorporated herein by refer- 
ence ) and a bit stream generated. This procedure is shown in detail in the flow- 
chart of Figure 5. representing a flow-chart of the P-frames coder, using atom 
based motion prediction. Dotted lines represent the coefficients and parameters of 
the atoms, the bold solid line represents a video frame and the light solid line 

10 represents bitstreams. The flow-chart should be understood as follow: The input 
frame are passed to the parameter estimation, which will modify the atom parame- 
ters stored in memory in order for them to describe correctly the input. The differ- 
ence between the new atom parameters and the one stored in memory is computed 
and the result is quantized and entropy coded. Those quantized difference are also 

15 de-quantized and added to the atom parameters previously stored in memory. This 
allows the reconstruction of the same atoms as the one available at the decoder. 
Those reconstructed atom parameters are then stored in memory, replacing the 
ones of the previous frame. They are also used to reconstruct the current frame 
using the generating mother function available at the encoder and decoder. The 

20 difference between this reconstruction and the original input is computed and en- 
coded using the frame coder. The number of updated atoms and their quantization 
can be fixed but can also be chosen in adaptive manner through rate and distortion 
constraints by means of a rate controller. For each motion predicted frame, the 
motion residuals (or texture) can be computed and encoded using any still image 

25 codec. It would then be sent in the bitstream in order to generate a scalable stream 
achieving lossy to lossless compression. 

Typical results of the motion prediction, with a coding cost of 170 Kbps in aver- 
age are shown in Figure 6, where predicted frames 20, 50 and 99 are represented. 
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It can be seen that the motion prediction stays accurate even after 100 frames, 
even in the absence of encoding of the motion residual. This is the major advan- 
tage of this technique compared to block based compensation. In order to show 
the advantage of the new prediction technique, we have done the same experiment 

5 with typical block matching compensation. The same sequence was encoded using 
adaptive block compensation with block sizes of 4x4 to 32x32. The encoder 
automatically selected the best block size according to Rate-Distortion optimisa- 
tion. By varying the block size it is possible to control the size of the compressed 
bitstream in order to match the result of the atom based motion prediction. Like in 

10 the case of the atom based prediction, the motion residual where not coded. The 
Motion prediction results of the block matching are shown in Figure 7. 

Objective comparison can be done using the PSNR measure. This measure is 
computed using the squared difference between the original and the reconstructed 

r \ 
255 2 



15 frame, i.e. = -lOlog 



\2 



X(/(x,y)-/ r (x,^)) : 



, where I(x,y) is the original 



frame and I r (x,y) is the reconstructed frame. The PSNR comparisons of the two 
prediction methods, for a bitstream of average size of 170 Kbps, are shown in 
Figure 8. The performance of the current atom based motion prediction is particu- 
larly good taken into consideration, that the states of the art block based motion 
20 prediction was used. Moreover, it can be seen that in the long term (more than 50 
frames) the prediction of the atom based method is more constant. In typical block 
matching applications, the block size is limited to 8x8, which causes poorer per- 
formances for the Block Matching. 



25 The number of predicted frames can be fixed but can also be chosen in adaptive 
manner through rate and distortion constraints by means of a rate controller. For 
instance, when abrupt changes occur (shots between scenes), the model is not an 
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accurate base anymore and the reference I-frame can simply be refreshed in order 
to compute a new model. Such a refresh mechanism can be monitored by tracking 
frame-to- frame distortion, which should stay rather constant for smooth updates of 
the initial model (Figure 8). 



